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Abstract. In this paper we develop a calculation code to account for the effects carried by fast rotation on the 
observed spectra of early-type stars. Stars are assumed to be in rigid rotation and the grid of plane-parallel 
model atmospheres used to represent the gravitational darkening are calculated by means of a non-LTE approach. 
Attention is paid on the relation between the apparent and parent non-rotating counterpart stellar fundamental 
parameters and apparent and true Fsini parameters as a function of the rotation rate stellar mass and 

inclination angle. It is shown that omission of gravitational darkening in the analysis of chemical abundances of 
CNO elements can produce systematic overestimation or underestimation, depending on the lines used, rotational 
rate and inclination angle. The proximity of Be stars to the critical rotation is re-discussed by correcting not only 
the Fsin i of 130 Be stars, but also their effective temperature and gravity to account for stellar rotationally 
induced geometrical distortion and for the concomitant gravitational darkening effect. We concluded that the 
increase of the Fsini estimate is accompanied by an even higher value of the stellar equatorial critical velocity, 
so that the most probable average rate of angular velocity of Be stars attains Q/Q c ~ 0.88. 

Key words. Stars: abundances - Stars: atmospheres - Stars: early-type - Stars: emission-line, Be - Stars: funda¬ 
mental parameters - Stars: rotation 


1. Introduction 


Fast rotation produces a polar flattening and an equato¬ 
rial stretching of stars, which in turn induce non-uniform 
surface gravity and tem peratur e distributions: the gravity 
darkening effect (GD, IZeipellll924 Sl. Owing to the lev¬ 
itation effect introduced by the rotation, rotating stars 
evolve as they had a lower mass: their core density is 
higher, the total bolometric luminosity produced is lower 
and they spend a longer time on the main sequence evo- 
lutionary phase than slowly rotating stars with the same 


mass 


lary phase than sfowiy rotating stars with the same 

llSackmannlll97(lt iBodenheimeill 1 97lt IClementlll97flt 


iMoss fc Smitlil ll982t IMevnet fc Maedejj 2000l) . The geo¬ 

metrical deformation of the star and the GD further in¬ 
fluence the apparent status of the star, which also de¬ 
pends on_wh£^erJt_Jias_aradiative_oi_^oiivectiveenve- 
lo pe | Luc vlll967t iTuominei] 19721 ISmith fc Worlevlll974t 


[a are H liook teonaT As predicted by models of stellar 

evolution tEnda l&L_Spfia| 1979t Ipinsonneault et al.lfl99lt 


iHeger fc Langeil l200(j iMevnet fc Maederl l200dh. fast 

tation is further expected to generate several instabili- 


ro- 
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ties that contribute to redistribute the internal angular 
momentum. They induc e thu s a turbulent diffusion in 
the stellar interior llZahnlll 992if . which in massive objects 
drives the CNO-cycled material from the core to the en¬ 
velope and ch anges the atmospheric che mical composition 
in these stars llMevnet fc Maedeill200f)l) . 

1.1. Rotation dependent spectra of early type stars 

Radiative and hydrostatic equilibrium conditions in atmo¬ 
spheres of rotating early-type stars impose in each point 
of the stellar surface the emitted bolometric flux be pro¬ 
portional to the local gravity F oc g^ s w ith 0 cs 1 i f T B ff A 
7000 K and g eS = g gra v - grot + g ra d iZeipe| l92j |Claretl 
1200011 . Although a non-local theory of radiative trans¬ 
fer in rotating stars does not reproduce von Zeipel’s re¬ 
lation, meridional circulation and turbule nce in surface 
layer s recover the validity of this law dSmith fc Worlevl 
[l974b . For B-type star effective temperatures, the radi¬ 
ation pressure g ra d can be neglected. Energy distribu¬ 
tions, photometric indices and spectral lines produced in 
gravity-darkened early type rotating stars were calculated 
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bv several authors. 

either for rigid 


1965, 

1966, 

1968a hlll974 

Hardoro & Strittmatterl 1968b, 

1972, 

1968? 

Hardoro & Scholz 197ll Gollins & Sonneborn 

1977 

1972 

Collins et al. 199lf Maeder & Peytremann 1970. 

or for differential rotators with conservative in- 


laws 


ternal rotational 
1986). AO to F5-type rigi d rotators were studied by 


Perez Hernandez et al ■ {ml). The effects of t he G D 


on 


the V sin i determinations were explored by Istoecklevl 
dl968h for rigid r otators and Zore c et al. (1988) for dif¬ 
ferential rotators. Istoecklevl lll968ll concluded that the re¬ 
duced contribution to the He I 4471 line broadening in the 
equatorial region may lead to underestimated rotational 
parameters in early-type stars, while Zorec et al. (1988a) 
obtained double valued Hsini in star m odels with polar 
hollows. Recently. iTownsend et al.l ( 2004 ) rediscussed the 
GD effect on the V sin i values of rigid early-type rotators 
and concluded that classic Rsinf determinations for BO 
to B9-type stars can be underestimated by 12 to 33% at 
u> = f2/f2 c = 0.95 rotational rates, if the Hel 4471 line is 
used and from 9 to 17%, if the Mgil 4481 line is studied. 


1.2. Aim of the present work 


WMgii (eq) > WMgii{ pole). Let us note that for the Hel 
4471 line, we actually mean the blend He I 4471+Hel 4470, 
where the second line is a forbidden component and has 
a somewhat different sensitivity to fundamental parame¬ 
ters than the permitted A4471 component. For the Hel 
4471 and Mgii 4481 lines also respond in a different way 
to the GD and because they overlap partially, the Hsinf 
determination demands that we proceed to fit both lines 
simultaneously. 

The aim of this paper is thus to take into account, as 
properly as possible, the changes of the stellar radii i? po i e 
and R eq and of the global bolometric luminosity with w 
in the calculation of the emitted radiation fluxes. We will 
also fit both Hel 4471 and Mgii 4481 lines in order to de¬ 
termine up to what degree the V sin i parameters of early- 
type fast rotators are underestimated if they are obtained 
assuming that rotational distortion and GD can be ne¬ 
glected. This is of particular interest for Be stars to know 
whet her they are near-critical rotators I Townsend et all 
In this paper, we will also discuss the effect of the 
rapid rigid rotation on the determination of the stellar 
effective temperature, surface gravity and chemical com¬ 
position. 


Apart from the obvious rotation dependent Doppler ef¬ 
fect, four factors, at least, still concur to produce the 
observed line broadening in gravity darkened fast rota¬ 
tors: a) changes with the surface non-uniform temperature 
of the continuous spectrum specific intensity; b) depen¬ 
dence of the intrinsic lines equivalent width with the local 
temperature and gravity; c) stretching of the isotempera¬ 
ture and isogravity regions, which make that in strips of 
constant radial velocity for rather large Doppler displace¬ 
ments there can be temperatures close to that of the undis¬ 
torted star; d) for high enough inclinations, the flattened 
stellar disc produce shortened constant Doppler displace¬ 
ment strips that tend to shallow the rotationally broad¬ 
ened line profile. A look on the Hel 4471 and Mgii 4481 
line equivalent width depe ndence with te mperature and 
gravity can be obtained in IPidelonl (Il982l) and Sect. n~3i 
These lines are currently use d to determine the V s in i pa¬ 
rameters of early-type stars dSlettebak et al.lfl975 ). From 
classic model atmospheres we obtain that at the He I 4471 
and Mg II 4481 line wavelengths the continuum specific 
intensity scales as I c ~ T“ ff , with a ~ 2.2 if T e g ^ 
7500 K. Since is it T eff (eq)/T eff (pole) = (g eq /g po ie ) 0 ' 25 , 
we infer that for uj = 0.95 the equatorial regions con¬ 
tribute less to the flux of the broadened line profile by a 
factor 7 c (eq)// c (pole) ^ 0.2 than in the pole. Line pro¬ 
files are then expected to be narrower than those cal¬ 
culated with models where GD is neglected. At high as¬ 
pect angles, effect d) compensates partially this deepen¬ 
ing. Depending on the temperature of the star, in fast ro¬ 
tators the change of the local T e g and log g with latitude 
produce the equivalent width of the Hel 4471 line be ei¬ 
ther W He i{e q) > W He i( pole), or W He i(e q) < W Hei ( pole), 
while the equivalent width of the Mg II 4481 line is always 


2. Model atmospheres 

Fast rotation flattens the star and produces non-uniform 
density and temperature distributions in its surface. To 
take into account the first order effects of this flattening on 
the stellar spectrum of B-type st ars, we adopted a me thod 
similar to the one described bv lGollins et al.l ( 199 11. but 
modified to introduce the changes detailed in Sect. 12.11 In 
this approach, that we applied in a computer code here¬ 
after named fastrot, the stellar photosphere is replaced 
by a mesh of plane-parallel model atmospheres depending 
each on the local temperature and surface gravity (Sect. 

IP) . 


2.1. Flattening and gravitational darkening 


We assume that the stars are rigid rotators without any 
surface latitudinal differential rotation component. We 
adopted the Roche approximation for the stellar surface 
equipotentials, i.e. in the gravity potential the dipole and 
higher order terms are neglected I Ta,ssonlll 978l) . In the ex¬ 
pression of the surface effective gravity we also neglected 
the latitude-dependent r adiation pressure term (g ra d = 0; 
iMaeder fc Mevnetll200Cll) . From the Roche equipotentials 
it follows that the ratio between the equatorial to the po¬ 
lar radii at the rotational rate u> = Fl/f 2 C ( fl c = critical 
angular velocity) is given by: 


-Re(^) 

R P (v) 

V 


1 + h \ 

1 


(i) 


where R c is the equatorial radius at critical rotation (w = 
!)■ 
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In what follows we use the fi/f l c to characterize the 
rotation rate instead of V e /V c where V e and V c are the 
actual and critical linear equatorial rotation velocities re¬ 
spectively. The main reason for this choice is that as we 
assume stars are rigid rotators, the angular velocity fl, 
which is the same over the whole star (surface and inte¬ 
rior) and independent of any other stellar quantity, it can 
be considered as an independent fundamental parameter. 
We note, however, that the stellar rotational distortion 
and the associated gravity darkening effect are controlled 
by the parameter 77 given in 0. Though f2/fl c , V e /V c 
and 77 are functions having the same germs 0 and 1 , it is 
f 2 /fl c < V e /V c < ?? in the open interval 0 < f 2 /f 2 c < 1 . 

In previous works it is often assumed tha t R„(uj) = R n 
wher e R n is the radius of the star at res t llCollins et alJ 
1 199 it iReinerslEooii iTownsend et al 1 12004 . In the present 
study we adopted instead the following interpolation ex¬ 
pr ession for . derive d from the calcul ation s performed 
bv iBodenheimeil l|l97llh IClementl (jl979ll and lZorec et, al.l 
<|l988tl : 

= l -P(M)t 

P(M) = 5.66 +^r, for M Z 2 

where M is the stellar mass in solar units and t = K/\W\ 
with K the rotational kinetic energy and W the stellar 
gravitational potential energy. A rough relation between r 
and 77 is given by r ~ [0.0072 + 0.008?; 1 / 2 ]?/ 1 / 2 . Although 
0 implies R c = 1.5 R 0 at critical rotation, as from m 
R p {to) < R a , depending on the B-type stellar mass, it is 
R c = 1.3R a to 1.4i? 0 , rather than R c = 1.5 R 0 as usually 
assumed in the literature. 

The total bolometric luminosity produced in the stellar 
core was assumed given by: 

L(uj,M,g 0 ) = L 0 (M,g 0 )f L (T,M) (3) 

where L a is the bolometric luminosity of an homologous 
non-rotating star with mass M and rest surface gravity 
g 0 . The expression for the factor /l was derived from 
Clement’s (1979) models for M <; 1.5: 



Relation © takes into account only the change of 
the total bolometric luminosity due to the mechanic ro- 
tationally induced c hanges of the st ellar core temper¬ 
ature and pressure |Saekma,nnl il97ctl . Whenever pos¬ 
sible, we have checked that the changes introduced by 
mixing phenomen a induced by the in s tabilities produced 
by the rotation l lEndal fc Sofial 1797(1 iMaeder fc Mevnefl 
ll996t iMevnet & Maeded 1200(1) do not introduce signifi¬ 
cant changes to the /l given by El 

In this work the GD law was written as: 

T e 4 ff (0,tu) = ' r (w,M,R o )g*(0,w) (5) 


where 9 is the colatitude angle; M is the stellar mass; R 0 
is the radius of the star at rest; g is the modulus of the 
gravity vector ~g = — V<f> with <I> being the Roche potential; 
/3 was considered = 1 for all local T e g (0, to) 7000 K. Fol¬ 
lower local T e ff, (3 was interpolated in Claret (1998). The 
factor 7 was calculated for each to as follows: 


L(u,M, g a ) 
crf s gPdS ' 


( 6 ) 


where a is the Stefan-Boltzmann constant and S is the 
total area of the surface Roche equipotential. From m 
it comes that the global effective temperature T e g = 
[L(to) / crS] 1 / 4 is smaller than T° s = [Lo/dTrcri? 2 ] 1 / 4 , the 
effective temperature of a non-rotating star with the same 
mass. On the other hand, relations and imply that 
TeffHpoie > T e ff. In the literature relation 0 is currently 
replaced by T 4 ff (0,w) = T 4 ff (0=O, u>)\g(0,uj)/g(6=0, w)]P, 
which hinders to relate the apparent stellar fundamental 
parameters with those the star would had if it were at 
rest. We note that the relations T e g(d=0 1 ui) = T° s and 
g(9 = 0,w) = g 0 , where T° s and g Q are the fundamen¬ 
tal parameters of the rotationless star, do not hold and 
that only the effective temperature and gravity averaged 
over the whole polar-on seen stellar hemisphere approach 
the effective temperature T° s and gravity g 0 roughly, i.e.: 
linq^o <T e g(6,u\i)> ~ T° s and lim^ 0 <g(0,u>\i)> ~ g Q . 
In all cases it is T e g(9= 0,co) > T° s and g(9=0, to) > g Q . 


/l = a{M) + [1 - a{M)\e~ b ^ T \ 

a(M) = 0.675 + 0.046M 1 / 2 \ (4) 

b(M) = 52.71 + 20.63M 1 / 2 J 

The quoted models (Bodenheimer 1971, Clement 1979, 

Zorec et al. 1988b) were calculated for stars with in¬ 
ternal conservative differential rotation laws and for en¬ 
ergy ratios r that reach the global secular stability limit 
r = A'/|kF| ~ 0.14 and even higher. In all these models 
the core is supposed to rotate as a rigid body. Moreover, 
these calculations show that the radii ratio R p (to)/R c and 
the function /l do not depe nd on the rotation law in the 
stellar envelope dZore lilii . Nevertheless, the interpola¬ 
tion relations 0 and 0 are for 0 < r = A7|IT| ^ 0.03, 
range of r values which are suited for rigid rotators, includ¬ 
ing the rigid critical rotation, where r = K/\W\ ^ 0.02 
whatever the stellar mass in the 2 ^ M ^ 60 range. 


2.2. Local plane-parallel models 

The spectrum of a rotating star is represented as the in¬ 
tensity of radiation emitted per unit wavelength interval 
per steradian in the direction towards the observer defined 
by the aspect angle i (inclination angle of the stellar rota- 
tion axis with respect to the line of s ight). It is given by 
( Collinsl fl965; IMaeder fc Peytreman S IlOZO ): 

r / 2 r n I,, I 

L x {i,u) = 2 &(j> R 2 (9)I x (p,oj)^-sm9d9 (7) 

Jo Jo cos $ 

where R(9) is the co-latitude 0-dependent distance be¬ 
tween the stellar center and the surface of the rotation 
ally-distorted star; p = p((f>,u}) = n. i (h is the unit 
vector normal to the stellar surface and i is the unit 
vector representing the direction of the line of sight); 
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cos 5(0, lu) = —h.f (r is the unit vector in the direction 
of R(6))] I\(fJ,,Lo) is the /^-dependent monochromatic spe¬ 
cific radiation intensity calculated for the local effective 
temperatures T e g(0,uj) and surface gravities g(9 1 u>). The 
integration of equation 0 was performed using a gauss- 
biquadrature of degree 96 in both angles 8 and </, which 
in principle is equivalent to a grid of 36481 surface el¬ 
ements of local plane-parallel model atmospheres. Each 
stellar surface point was characterized by its specific ra¬ 
dial velocity in the observer’s direction in order to take 
into account the corresponding Doppler shift in the spec¬ 
tral line. 


The models we used to evaluate each local spe¬ 
cific intensity I\{p,u>) were computed in two consecutive 
steps. To account in the most effective way for the line¬ 
blanketing, the temperature structure of the atmospheres 
have been computed as in Kurucz fc Cast elli (2000) using 
the ATLAS9 computer code (lKurucd|l993l) . Non-LTE level 
populations were then comp uted for each of the a toms 
we considered using tlusty llHubenv fc Tanzlll 99ft) and 
keeping fixed the temperature and density distributions 
obtained with ATLAS9. 


Excepted for C II, the atomic models we used in 
this work were downloaded from TLUSTY’s homepage 
(http://tlusty.gsfc.nasa.gov) maintained by I. Hubeny and 
T. Lanz. Table^lists the ions that were introduced in our 
computations. C n was trea ted w ith use of the MODION 
IDL package developed by IVarosi et alJ lll995l) and by 
adopting the atomic data (oscillator strengths, energy lev¬ 
els and photoioniza tion cross-sections ) selected from the 
topbase database dCunto et all 19931) . It reproduces the 
results obtained bv ISigurt lllOOfit) . 

In this way and for each spectral region studied in the 
present work, the specific intensity grids were computed 
for effective temperatures and surface gravities ranging re¬ 
spectively from 15000 K to 27000 K and from 3.0 to 4.5. 
For T° s < 15000 K and T° s > 27000 K ( T° s concern rota¬ 
tionless model atmospheres) we used LTE calcu lations and 
the O STAR.2002 model atmospheres grid dhanz fc Hubenvl 
12003 ). 


3. Comparison of rotating with non-rotating 
models 

3.1. Adopted procedure 

Hydrogen and helium line profiles, more particularly, the 
Hy, A4388 Hei and A4471 Hei lines, are often used to 
derive the effective temperature, surface gravity and pro¬ 
jected rotation velocity of B type stars. Among the rea¬ 
sons of their use, we have: 1) the line profiles and equiva¬ 
lent widths are very sensitive to effective temperature and 
surface gravity; 2) the blue helium lines are known to be 
somewhat less sensitive to LTE departures; 3) the broad¬ 
ening mechanisms of thei r line profiles have been well 
known for a long time (e.g.: B^rn^n^et ; all19691 iLeckron d 


Table 1 . Atomic models used 


Atom 

Ion 

Description 

Hydrogen 

H I 

8 levels + 1 superlevel 


H II 

1 level 

Helium 

He i 

24 levels 


He ii 

20 levels 


He hi 

1 level 

Carbon 

C ii 

53 levels all individual levels 


C hi 

12 levels 


C iv 

9 levels + 4 superlevels 


C v 

1 level 

Nitrogen 

N i 

13 levels 


N ii 

35 levels + 14 superlevels 


N hi 

11 levels 


N iv 

1 level 

Oxygen 

O i 

14 levels + 8 superlevels 


O ii 

36 levels + 12 superlevels 


O III 

9 levels 


O IV 

1 level 

Magnesium 

Mg ii 

21 levels + 4 superlevels 


Mg hi 

1 level 


Table 2. Sets of parameters used in model grids 


GRIDA 

5000 K < 

rjiO 

1 eff 

< 

40000 K 

step = 

1000 K 

2.5 < 

log So 

< 

5.0 

step = 

0.5 

0 km s -1 < 

Esin i 

< 

500 km s _1 

step = 

50 km s _1 

GRIDB 

12000 K < 

rjiO 

eff 

< 

27000 K 

step = 

1000 K 

3.2 < 

logfito 

< 

4.4 

step = 

0.2 

0° < 

i 

< 

90° 

step = 

5° 

0.50 < 

n/n c 

< 

0.99 

step = 

0.05 


Two grids of spectra ranging from 4250 A to 4500 
A were computed. A first grid, GRIDA, was obtained 
using classical plane-parallel model atmospheres. A sec¬ 
ond grid, GRIDB, was calculated using fastrot (Sect.0 
where stellar flattening and gravity darkening are taken 
into account. The sets of fundamental parameters consid¬ 
ered in these grids are given in Tabled T° s and g Q refer to 
the fundamental parameters of the rotationless or parent 
non-rotating counterpart of the star f Sect. mi) . Regarding 
the notation, it is worth noting that for GRIDA it holds 
T e s = T° s and log g = log-¬ 
in what follows we will determine the effective temper¬ 
ature, surface gravity and the projected rotational veloc¬ 
ity of B-type stars using only the Hy, Hei4388, Hel4471 
and Mgli4481 lines. To this end we will use a fitting pro¬ 
cedure based on a least-squares method that applies the 
MINUIT minimization package available at CERN. In or¬ 
der to know the confidence degree of this method, we 
fitted the spectra of the GRIDA reference lybrary using 
synthetic spectra obtained from the same NLTE plane- 
parallel model atmospheres like the GRIDA spectra. The 
free, or fit parameters used to this purpose were: the ef- 
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fitted log g fitted V sin /( km s* 1 ) 


Fig. 1. Test of the adopted procedure: comparison of T e s , log g and Fsini parameters obtained by the fitting procedure 
with the model parameters of GRIDA 


Table 3. Spectral regions considered for the fits. 


Spectral line 

CASE 

A B 

AA (A) 

Hy 

X 

X 

4300 - 4375 

A4388 He i 

X 

X 

4375 - 4400 

A4471 He i 


X 

4460 - 4475 

A4471 He i + A4481 Mgn 

X 


4460 - 4490 


fective temperature, surface gravity, the projected rota¬ 
tional velocity and the averaged flux ratio of the fitted to 
the reference spectrum. The \ 2 deviation parameter was 
computed only for the selected spectral domains defined 
in Table[o]for CASE A. The comparison shown in Fig. □of 
parameters obtained by the minimization procedure with 
the model GRIDA reference parameters, reveals that the 
agreement is quite satisfactory. However, this agreement 
is the best for the projected rotational velocity, whose ac¬ 
curacy is generally better than 4%. The derived effective 
temperature approaches the most to the reference values 
when T e ff ~ 20000 K, but for other temperatures devia¬ 
tions never exceed 500 K. The errors committed on logg 
are always lower than 0.05 dex whatever the effective tem¬ 
perature. 

To study the effects of fast rotation on the determi¬ 
nation of stellar fundamental parameters, we applied the 
same procedure than above to the GRIDB spectra, where 
the model parameters are {T° ff , log g Q , i,2l/H c } or equiv¬ 
alently {T^,logg 0 ,Fsini,f2/fI c }. Since the fit of spectra 
produced by rotating stars can be done in two ways, we 
distinguish two sets of parameters. The reference model 
parameters proper to calculate fastrot rotationally mod¬ 
ified synthetic spectra, as well as those obtained by fitting 
with them the rotationally distorted spectra, are hereafter 
called parent non-rotating counterpart (pnrc) fundamen¬ 
tal parameters and will be noted with a superscript “ ° ” 
or a subscript “ 0 ” (i.e. T° s and logg 0 ), while Fsini will 
represent the true projected rotation velocity of the star. 
Conversely, we call apparent fundamental parameters (i.e. 
T e a g p ’, log <?app. and V sin i app .), those derived by fitting the 


GRIDB spectral lines with synthetic spectra issued from 
classical plane-parallel model atmospheres. 

The fit of the Hy, Hei4388 and Hel4471 lines was done 
either including the Mgn4481 line or excluding it in order 
to test the procedure’s sensitivity to a change of the fitting 
criteria; this identifies respectively CASE A or CASE B 
in Table 0 


3.2. Reference parameters 


In the following sections, calculations will be done for sev¬ 
eral effective temperatures and surface gravities. To have 
a rough look on the spectral types they may concern as 
well as the relations that exist between the fundamental 
parameters suited for non-rotating stars and the apparent 
ones - which are dependent on the stellar geometrical dis¬ 
tortions and the concomitant GD - we list in Table 01 the 
reference parameters we used in the present work. In this 
table are also given the polar effective temperature and 
the polar radius of the distorted stars for several values of 
the angular velocity ratio f l/fl c . ^ 

The MK spectral types adopted I Gray fc Corballvl 
imi are for the effective temperatures and logg 0 = 4.0 
of the parent non-rotating stars. The stellar masses and 
radii for S2/f2 c = 0 used to calculate the critical velocity 
V c (c.f Sect. 12 . Ill were deri ved from stella r evo lutionary 
tracks of non-rotating stars llSchaller et al.lll99 2b 


4. Apparent fundamental parameters 

Following the procedure described in Section cm we ob¬ 
tained a grid of pnrc (i.e.: parameters belonging to the par¬ 
ent non-rotating stellar counter-part) and apparent fun¬ 
damental parameters for different values of the angular 
velocity and inclination angle of the star. Since the aim 
of the present computations is to provide a way to cor¬ 
rect the observed fundamental parameters from gravita¬ 
tional darkening effects for a large parameter-space, in the 
present work we have chosen to present the effects due 
to fast rotation and the corresponding corrections, as a 
function of T a | p -, log 5 app . and Fsini app . (e.g. Fig. 0, 
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Fig. 2. Apparent effective temperature and surface gravity as a function of Vsini for different pnrc (T e ° ff ,logg 0 ) values 
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Table 4. Polar temperature (T po i ar ), polar radii (i? p ) 
and critical velocities (E c ) as a function of f2/f2 c and for 

log 9o~ 4. 


MK Type 

rriO 

eff 

K 

Mass 

M© 

P/Q c 

-tpoiar 

K 

R P 

R© 

Ec 

km/s 

B1 V 

24620 

10.53 

0.00 

24620 

5.32 

525 




0.50 

24932 

5.22 





0.80 

25887 

5.10 





0.99 

27593 

4.88 


B2 V 

19500 

6.87 

0.00 

19500 

4.30 

472 




0.50 

19736 

4.21 





0.80 

20471 

4.11 





0.99 

21773 

3.94 


B5 V 

14000 

4.15 

0.00 

14000 

3.34 

418 




0.50 

14176 

3.27 





0.80 

14708 

3.19 





0.99 

15609 

3.04 


B9 V 

10720 

2.86 

0.00 

10720 

2.77 

383 




0.50 

10870 

2.71 





0.80 

11300 

2.63 





0.99 

12007 

2.50 



Table 5. Effects of fast rotation in equator-on B 
type stars. The effective temperature (RT e g) and sur¬ 
face gravity (Rlog g) underestimation percentage are 
given against f 2 /f 2 c . 


n/fic 

RTeff 

Rlog g 

0.80 

9 % 

8 % 

0.90 

12 % 

10 % 

0.99 

17 % 

20 % 


to make a direct link with the actually measured quanti¬ 
ties, rather than against the pnrc fundamental parameters, 
which need to be determined from the sought corrections. 
As the REFs are strongly dependent on the aspect angle *, 
we present them as a function of the observed Esin i app . 

4.1. Effective temperature 

The typical effects of stellar flattening and GD on effective 
temperature are shown in Fig.|2]for different combinations 
of T° s and log<? 0 , while the corrections A T e g to add to the 
apparent effective temperature to get the pnrc effective 
temperature at different apparent f 2 /f 2 c , apparent Esin*, 
and apparent log g regimes are plotted in Fig. 0 (upper 
panel). As expected, the T^ p ' decreases with increasing 
Esin* app .. Fast rotation produces therefore a lowering of 
the apparent effective temperature. Concerning the ampli¬ 
tude of this effect and its general behaviour, the following 
remarks can be made: 

i) REFs on T a g P ' become significant for fi/f2 C > 0.6. 

They are stronger and vary more rapidly with 

Esin* app . the larger the values of f l. 


**) When fi/increases, REFs on T a g P ’ become more 
sensitive to surface gravity in the sense that they are 
larger the lower is the gravity. 

Hi) The relative lowering effect on the apparent T e g does 
not depend much on its pnrc value. In Table 0 this 
effect is given in percentages for different angular ve¬ 
locities in stars seen equator-on, where effects are the 
stronger. 

iv) When studying dwarf B-type stars, REFs on T a g P ' can 
be ignored safely at Esini app . lower than 200 km s _1 
and at Esin* app . ^ 100 km s -1 for giant B-type stars. 

4.2. Surface gravity 

Corrections Alog g to add to the apparent surface grav¬ 
ity to get the pnrc surface gravity at different ap¬ 

parent Esin* and apparent T e g regimes are plotted in 
Fig- El (lower panel). As shown in Fig. El stellar flattening 
and gravitational darkening produce a lowering of the ap¬ 
parent surface gravity. The amplitude of this REF varies 
obeying the following rules: 

*) REFs on loggapp. become significant for >0.6. 

They are stronger and vary more rapidly with 
Esin* app. at high values. 

**) When ^ 0.85, REFs on logg ap p. are weaker in 
stars with T° s 20000 K than in those with T° e 5s 
20000 K. On the contrary, for H/fl c <; 0.85, the appar¬ 
ent log g are more sensitive to REFs in stars with T° s 
20000 K than in those with T° s ^ 20000 K. This 
change in the log g sensitivity to REFs, according to 
T° s and is due to hydrogen lines, which in the 

present approach are the main log g indicators. In hot 
stars, the intensity of hydrogen lines is the lowest, but 
it increases fast with the equatorial local T e g decrease 
produced by increasing values of fi/fl 2 C . 

Hi) The lowering of log Sapp, does not depend much on the 
pnrc log g value. It is given in percentages in Table 
El for different angular velocities in stars seen equator- 
on where effects are the stronger. In the more distorted 
cases, the change in logg a pp. represents two luminosity 
classes. 

iv) In B-type stars with T^ p ' 5 ) 18000 K and Esin* app . 
^ 100 km s -1 , REFs on log g can be ignored. 

4.3. Projected rotation velocity 

The REFs on Esin* ap p. become significant for f2/f2 c <; 
0.70. The spectral lines become progressively less broad¬ 
ened than expected from models of stars without GD, 
due to a lower contribution to the line flux coming from 
the gravity darkened equatorial regions towards the stellar 
limb. The magnitude of this effect depends on the sensitiv¬ 
ity of the line studied to the local changes of T e g and log g , 
as shown in Fig.EJfor the Hei 4471 and for the Mgn 4481 
lines. 

At a given f2/f2 c 0.7 and T° s 18000 K, the 
broadening of the Hei 4471 line tends to saturate for 








A log (g ) at £i/£i c =0.80 A log (g ) at £i/£i c =0.99 
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Fig. 4. Equivalent widths of the He I 4481 and Mgn 4481 
spectral lines against effective temperature for different 
log g values. Computations were made assuming plane- 
parallel model atmospheres. 



Fig. 5. Apparent Fsinz derived using the Fourier 
Transform technique against the true Fsin* considering 
individually the Hei 4471 line (panel a) and Mgli 4481 
spectral line (panel b). Computations are made assuming 

n/n c = o.99. 


increasing values of Fsiniapp.. In stars cooler than T° s 
~ 18000 K, the saturation of this broadening is favored 
also by the enhanced contribution to the line flux from 
the polar regions, where the local T e g is increased by the 
CD. Fig. 0 shows the apparent V sin zfft against true 
Fsinz for logg 0 = 4.0, fi/S7 c = 0.99 and for different 
values of T° s when only the He I 4471 line (Fig. 0 a) is 
considered. The apparent Fsinz parameter is obtained in 
this figure using the Fourier transform method applied 
to line profiles calculated with FASTROT. In the present 
work, the Fourier method is meant to represent the clas¬ 
sical techniques of determining rotational parameters. We 
see that in a star having a T° s = 14000 K, a Fsini app . 
~ 250 km s' 1 value implies that the underestimation 
of the projected rotational velocity is (5Fsini= (Vsinz)- 
(Fsinz app .)— 140 km s -1 , if it rotates at H/H c =0.99. 
The underestimation SV sin i drops fast as soon as H/H c 
< 0.99, while the magnitude of the Vsini app saturation 
depends on the inclination angle, angular velocity and also 
on the effective temperature (e.g. Fig. 000 and 0). 


Same conclusions can be drawn for the Mgn 4481 line. 
However, in this case the above mentioned line broaden¬ 
ing underestimation is less pronounced (Fig. 0b). This is 
due to the fact that the intensity of the Mgn 4481 line 
is smaller in the polar regions than in the equatorial, so 
that even the continuum intensity in the line wavelengths 
is stronger at the pole than in the equator. The differ¬ 
entiated response of the continuum and the intrinsic line 
intensity to T e g- and log g compensate to each other some¬ 
what, producing hence a smaller 5Fsin i. 

The relation between the apparent and the true 
V sin i derived from the Mgli 4481 is shown in Fig.0 We 
note in this figure that saturation effects are seen in the 
He I 4471 line for low T° s values and in the Mgn 4481 
line when T° s is high. In the first case, the effect is due to 
a significant decrease of the equivalent width of the line 
for low effective temperatures (i.e. at the stellar equator). 
In the second case, the contrast of the continuum specific 
intensity between the pole and the equator cannot be com¬ 
pensated by the modest increase of the equivalent width 
of the line at lower equatorial temperatures (cf. Fig. 0. 
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V sin / ( km/s ) 


Fig. 6. Comparison between true and apparent Vsinz in 
fast rotating stars assuming logg app .= 4. The fitted re¬ 
gions correspond to CASE A of Table 0 


As can be seen in Fig. El where Vsinz app . was de¬ 
rived following the procedure described in Sect. 3.1, the 
simultaneous fitting of various spectral lines combined 
with the simultaneous determination of T^g P ', logg app . 
and Vsinz app limits somewhat the saturation effect and 
minimizes the underestimation of the projected rotation 
velocity at high inclinations and angular velocity rates. 
Fig. 0 shows the relation between the apparent and true 
Vsinz for f2/f2 c = 0.80, fi/fi l c = 0.99, logg 0 = 4.0 and for 
different T° s values when both He I 4471 and Mgn 4481 
lines are taken into account for the Vsinz determination. 
Considering the case f2/f2 c = 0.99 (i.e. V e /V c — 0.97), we 
see that for whatever pnrc T e ff, the apparent Vsinz under¬ 
estimates the true rotation by SV sinz ^ 50 km s ^ 1 as far 
as the apparent Vsinz 350 km s -1 . Underestimations 
may then increase to ('Vsinz ~ 80 km s _1 for apparent 
Vsinz 350 km s -1 . It is worth noting the peculiar be¬ 
havior of this relation for the T° s = 15000 K. The en¬ 
hanced contribution of the Mgil4481 line to the blend 
Hel4471+Mgll4481 at high inclinations, makes that the 
Vsin z is likely determined by the behavior of the Mg II4481 
which, as seen in Fig. El leads to a smaller SV sinz than the 



. 

0.5 0.6 0.7 0.8 0.9 1.0 

V e sin //V e 

Fig. 7. Apparent and true Vsin z values at T° ff = 20000 K, 
l°gffo= 4.0 and different values of sin z (full-line curves) 
and of H/f 2 c (dotted curves). 


He I line alone. It is then crucial that for stars with T° s <, 
15000 K both Hel4471 and Mgil4481 are considered to¬ 
gether in the Vsinz calculation. On the contrary, for T° s 
*5 20000 K, the addition of A4481 Mgn in the fit does not 
increase significantly the accuracy of the Vsinz measure¬ 
ments. Fig. El also shows that for high enough H/H c ratios 
the underestimation SV sinz is a function of T° ff . 

As the use of the Vsinz parameter at a given H/H c 
hides the effect of the aspect angle on the measured Vsinz 
we plotted in Fig. Q for a pnrc B2 V spectral type, the 
behaviour of the apparent V sin z against its true value at 
different fixed inclination angles (sinz). We note that the 
underestimation at higher Vsinz is also very sensitive to 
the sinz. 

The incidence of log g Q on the underestimation SV sin z 
is further shown in Fig. El In this figure the relation be¬ 
tween i5Usin z against the true V sin z is given for differ¬ 
ent pnrc effective temperatures and surface gravities. It 
is seen that as long as Vsinz <j 300 km s -1 , the JUsinz 
is almost independent of gravity. As soon as Vsinz 300 
km s~\ not only do the SV sin z differences depend strongly 
on log go, but depending on the effective temperature it 
shows two different behaviors. For temperatures T° s ^ 
20000 K, <5Vsinz begins to decrease, while for T° s <; 20000 
K it increases even faster. In Fig. El the curve correspond¬ 
ing to log go = 3.4 is badly determined for high Vsinz 
values, because the gravity decreases rapidly to very low 
values at the equator. 
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Fig. 8. Vsini lowering (<5Ksini) as a function of the true 
Vsini for different true effective temperature and true 
log g values. 

5. Effects of fast rotation in other spectral regions 

It is very important to know whether plane parallel model 
atmospheres calculated for the apparent fundamental pa¬ 
rameters can describe accurately enough the whole ob¬ 
served spectrum of a rotating star. This matters not only 
for the study of the stellar chemical composition, but 
also for the energy distribution representation, which en¬ 
ters frequently the fundamental parameter determination. 
This is important, in particular, for fast rotating early 
type stars, whose evolution is faster than that of less mas¬ 
sive stars, so that they are found in different evolutionary 
stages and can be used to test models of stellar evolution. 

5.1. Chemical composition 

In order to verify the reliability of the use of computed 
stellar spectra that ignore GD, we studied wavelength do¬ 
mains that contain several helium and hydrogen lines, as 
well as CNO spectral lines. We compared apparent spec¬ 


tra, i.e. those computed with a plane-parallel model atmo¬ 
sphere using the apparent fundamental parameters given 
in Tabled with pnrc spectra (i.e. those calculated using 
fastrot for T° s =20000 K, logg 0 =4) and several incli¬ 
nation angles f2/Q c = 0.99. This comparison for hydrogen 
and helium lines is shown in Fig.d The same type of com¬ 
parison for CNO elements is shown in Fig. m It appears 
then that in general REFs can be ignored for stars with 
Vsini^l50, which correspond either to slow rotators or 
pole-on fast-rotators. However, in fast rotating stars seen 
at high inclination angles, ie: large Vsini values, discrep¬ 
ancies between both series of models can be significant. 
Thus, the chemical composition of objects seen at high in¬ 
clinations, which are also the most frequently found, must 
be studied using models adapted for rotating stars. 

In panel (a) of Fig.[^]are shown the helium and hydro¬ 
gen lines used to determine the apparent T e s and log g , 
as noted in Sect. rm The spectral regions where the fit 
was actually performed are marked on top of spectra with 
horizontal thick lines. In general, these same plane-parallel 
models calculated and apparent fundamental parameters 
are representing quite satisfactorily the helium and hydro¬ 
gen lines in other spectral regions (see panel (b) in Fig. |9j) . 

The behaviour of the CNO spectral lines as function of 
the inclination angle and angular velocity is however more 
complex, as shown in Fig.EB In Fig. mi we reported the 
overabundance of oxygen and nitrogen derived for models 
corresponding to an apparent B2 IV star (T2 j P = 19000 K 
and log(/app.=3.5) but for different values of and 

two different Vsini values (150 km s -1 and 200 km s^ 1 ). 
The adopted fundamental parameters are listed in Tabled 
This kind of modeling shows that in early B-type stars 
C II is expected to be underestimated when use is made 
of plane-parallel model atmosphere for fast rotators seen 
equator-on. Nevertheless, these models can still be used 
when stars are seen at low inclinations. The abundance 
of nitrogen, as derived from various N n lines, is overes¬ 
timated by a plane-parallel modeling of fast rotators seen 
at intermediate inclination angles. In fast rotating B-type 
stars cooler than 25000 K, oxygen is always overestimated 
by plane-parallel models when fits of On spectral lines are 
done, while it is expected to be somewhat underestimated 
if one uses the A7772 Oi triplet. 

More generally, spectral lines preferably formed at the 
stellar poles (N II, O II ...) tend to appear stronger when 
the effects of fast rotation are taken into account (see f.e. 
Oil lines in panel (a) of Fig. 1101) . On the contrary, when 
their privileged forming region is near the equator, they 
often appear weaker (see f.e. Oi line in panel (b) of Fig. 
HU but the REFs are also generally smaller. Therefore, 
REFs on chemical abundance determinations strongly de¬ 
pend on the studied transition, on Vsinz and on the effec¬ 
tive temperature of stars. The content of CNO elements in 
atmospheres of early-type fast rotators, in particular Be 
stars, will be discussed in a forthcoming paper. 
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Fig. 11. Abundance overestimation in dex due to 
stellar flattening and gravitational darkening observed 
for different values of Q/Q c and Vsini app . equal to 
150 km s -1 (open circles) and 200 km s -1 (filled circles). 
Adopted fundamental parameters are listed in Table ||3 
Spectral lines used in the abundance estimate are noted 
in each Figure. 




1974]; Hardorp & Strittmatter 1968a. 1972, 

1968a 

Hardorp & Scholz 

197ll Collins & Sonneborn 

Il977 

Collins et, al. 199l| 

Maeder & Pevtremann 11970, 

1972 


Cq^u^^^mjTk|j_fl85[|Zore^ fl08(f|p^re7Hern5rKiey et aiJ 

19991 iTownsend et al.ll20041. To complete the presen¬ 


tation of rotational effects and because the spatial 
UV spectra l reg ion was less e xplo red in the literature 
jKodaira & Hoekstralll979t iLlorente de Andres & Du rail 
Il97f)l l we show in Fig.fT 2 l the comparison between spatial 
UV energy distribution computed with FASTROT using 
true and apparent parameters for U/U c = 0.99 and 
different inclination angles. The values of parameters 
used in each calculation are given in Table Q Sometimes 
the spatial region of stellar energy distributions are used 
t o compare appare nt diameters or to infer distances 
( Fremat et alJ 1200 2). According to estimates shown in 
Fig. El differences up to 70% are expected in these 
quantities if they concern early type fast rotators. 


Table 6. Apparent (T e a | p ', logg app ., Vsini app .), pure (T° ff , 
log go) parameters and true Usini used to study the ef¬ 
fects of stellar flattening on the determination of CNO 
abundances IFig.ITTl). 


cy-iapp. 

1 efi 

(K) 

l°g Sapp. 

Vsini app . 
(km s -1 ) 

Q/R c 

rjiO 

eff 

(K) 

log So 

V sin i 
(km s -1 ) 

(°) 

19000. 

3.5 

150 

0.80 

19433 

3.64 

158 

38 

19000. 

3.5 

150 

0.85 

19709 

3.62 

160 

34 

19000. 

3.5 

150 

0.90 

19662 

3.65 

161 

30 

19000. 

3.5 

150 

0.95 

19537 

3.66 

164 

27 

19000. 

3.5 

150 

0.99 

20480 

3.74 

170 

24 

19000. 

3.5 

200 

0.80 

19529 

3.66 

207 

50 

19000. 

3.5 

200 

0.85 

19680 

3.69 

210 

45 

19000. 

3.5 

200 

0.90 

20070 

3.71 

212 

38 

19000. 

3.5 

200 

0.95 

20424 

3.75 

218 

35 

19000. 

3.5 

200 

0.99 

20818 

3.75 

227 

32 


Table 7. True and apparent fundamental parameters used 
in Fig. El and Elat U/fl c =0.99 . 


i 

rpO 

eff 

log So 

Vsini 

j-iapp. 

J eff 

log Sapp. 

Vsin i app. 

0 

20000 

4.0 

0 

18675 

3.83 

0 

30 

20000 

4.0 

231 

18462 

3.75 

203 

90 

20000 

4.0 

463 

16249 

3.28 

382 


5.2. Energy distribution 

Since REFs on the stellar continuum are directly pro¬ 
portional to the shape and apparent size of the observed 
stellar disk, magnitudes are very sensitive to the stellar 
flattening. The GD displayed by the stellar hemisphere 
facing the oberver enhances even more these effects. REFs 
affect therefore the slope and the absolute level of the 
emitted fluxes. REFs on the stellar continuum in the vis¬ 



Wavelength ( A) 

Fig. 12. Comparison between spectra computed using 
pnre (full line) and apparent parameters (crosses). Pnrc 
spectra were computed for T° ff = 20000 K, logg 0 = 4.0 and 
fl/n c = 0.99. Three different inclination angles were con¬ 
sidered. Apparent values for each case are given in Table0 


ible spectral range ha ve already been widely discussed in 
many previous works JCollmdll 96.4 Il965l Il flfid Il968alhl. 
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Fig. 13. Comparison between Hei A4471 line profiles com¬ 
puted with BRUCE (full line, Townsend et al. 2004) and 
with FASTROT (crosses). Adopted input parameters in 
FASTROT are: T° s = 20000 K, log ffo = 4.0, U/U c = 0.99 
and Usin i =463 km s -1 (ie: polar temperature = 22300 K 
and polar radius = 4.03 R©). Adopted input parameters 
in BRUCE are: polar temperature = 22300 K, polar ra¬ 
dius = 4.03 RO, inclination angle = 90° and equatorial 
velocity = 463 km s _1 . 


6. Discussion 


In the present approach of rotational effects on spec¬ 
tra emitted by early-type stars we neglected second or¬ 
der effects due to radiative fluxes related to the latitu¬ 
dinal gradient of the local effec tive temperatu re. They 


nperature. 

dom, i 


might be not entirely negligible iHadravali I uwa i . m par¬ 
ticular for very fast rotators, where the equatorial gravity 
lim^^i g(e q) —> 0. It is expected that their effect may 
reduce somewhat the GD. 


6.1. Comparison with other recent calculations 

To test our code, we compare in Fig. El the Hei 
A4471 line profile com puted with FASTROT and bruce 
iTownsend et al.l2004lf . The same model atmosphere grids 
were used assuming near critical rotation f2/f2 c = 0.99 
(V e /V c = 0.97) and inclination i = 90°. Relations 0 
and m were combined to find what pnrc effective tem¬ 
perature is implied by the T po i ar value we used as entry 
parameter in bruce. This procedure leads to an almost 
identical gravitational darkening and to identical line pro¬ 
files as shown in Fig. d We made several comparisons of 
FASTROT and BRUCE results of the type shown in Fig. ITTH 
Although our integration algorithms is different and out¬ 
numbers the equivalent elementary surface elements those 
of BRUCE, quite similar results are obtained. The only 
difference between both approaches relies on the adopted 
fundamental parameters. In FASTROT a unique set of pnrc 
T° s and mass M is used for the whole series of U values, 
while BRUCE does allow to specify the temperature distri¬ 


bution in 2 possible ways: by fixing_T afl j al _or by fixin g th e 
bolometric luminosity (see eq. 3 in iTownsend et al . 2j)0|. 


6.2. Implications for the study of Be stars 


It has been pointed out by IStoecklevl (Il968t) and 
ITownsend et alJ ll2004tl that the FWHM of the A4471 He I 
line can converge to a constant value at high V sin i when 
the stellar angular velocity approaches the critical rota¬ 
tion, which makes that the Usinz values of fast rota¬ 
tors, in particular those of Be stars, can be underesti¬ 
mated. Nevertheless, the implication on the measurements 
of V sin i in actual early-type stars by taking into account 
entirely the blend Hei4470+Hei4471-|-Mgil4481, has still 
not been entirely explored. Our comparisons between true 
and apparent values obtained using the fitting procedure 
(CASE A Table 0, which mimic real cases of V sin i deter¬ 
mination, show that there is a progressive lowering of the 
Usin * for increasing rotational rates. We also see that this 
V sin i underestimation does not imply necessarily a satu¬ 
ration of the estimated projected rotation velocity, even 
there is critical rotation (see Fig. 0 . 


It is interesting then to revisit the most probable value 
of the fl /flr rate of Be stars if the REFs are taken into 


Chauville et al . fcomh have studied a sample of 


account. 

116 Be stars and determined their Usini app . parameter 
without taking into account the REFs. They used, how¬ 
ever, Stoeckley and Mihalas’ (1973) calculations of rota- 
tionally broadened lines, where the wavelength dependent 
limb-darkening within the line profiles is considered and 
which ensures that there is no systematic underestima- 
tio n of rotational p a ramet e rs, effect that was dis cussed 
by ICollins fc Truaxl lll99fil) . IChauville et al.l ll200llf have 
thus concluded that on average Be stars rotate at U/U c ~ 
0.8. We have studied these stars again, and some others, 
and redetermined carefully their apparent fundamental pa¬ 
rameters by fitting their spectra with model atmospheres 
without r otation. We adopted t he apparent Usinz deter¬ 
mined bv I Chauville et all d200lh in order to avoid unde¬ 
terminations of these quantity produced by the Fourier 
method, which uses a constant limb-darkening coefficient 
in the frequencies over the line profile. Using our GRIDB 
of synthetic spectra, we have determined the pnrc fun¬ 
damental parameters of the program objects, in partic¬ 
ular their Usinz. These parameters were calculated for 
an average fi/f2 c rate, the same for all stars and whose 
value was iterated. The iteration started first by assuming 
f2/f2 c = 0.8. At each iteration step we corrected the appar¬ 
ent fundamental parameters in order to obtain the pnrc 
ones and redetermined, hence, the ratios Usinz/U, of each 
stars and studied its distribution. Each time we derived 
from the distribution the average of fl/Q c and its mode, 
i.e. the most frequent or probable value of S2/U c . These 
quantities were estimated in two ways: a) using the same 
method as IChauville et al ■ (l200ll) : b) determining the av¬ 
erage and the mode of V/V c by correcting the distribution 
of Usinz/V), for the sinz effect. Since the mode of a given 
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distribution depends on at least its first four moments, 
where errors increase with the order of the moment, we 
iterated the average We also started the iteration 

assuming all stars rotate at f2/f2 c = 1.0. We note that each 
apparent effective temperature, gravity and projected ro¬ 
tational velocity was considered with its 2 a error. The in¬ 
terval (X — 2ax, X + 2i tx) corresponding to each apparent 
X parameter was divided into eight parts, so that we had 
9 entries for each independent parameter. This produced 
9 3 determinations of each corresponding pure X parame¬ 
ter or true Tasini. Since each distribution of the obtained 
pnre parameters is not symmetrical, we adopted its mode 
to represent the corresponding most probable solution of 
the sought pnre parameter. 

Both iterations ended at the same values of the average 
and the mode: 

= 0.85±g;g| (fi/fi c ) m «fc = 0.88ig;jg (8) 

Table E summarizes the characteristics of iterations. 
Moments // n are for V/V c distributions. The program stars 
and their fundamental parameters are given in Table El In 
this table are given the apparent effective temperature, 
surface gravity and rotational parameter with their re¬ 
spective la errors. The modes of the corresponding pnre 
parameter and of the corrected I^sin* are also listed, as 
well as the critical equatorial velocity, the estimated in¬ 
clination angle and the ler dispersion of the respective 9 3 
determinations of each fundamental parameter. 

From the rates displayed in (0 we see that critical ro¬ 
tation among Be stars is reachable only by a fraction of 
objects between 2 and 3cr interval beyond the ( Q/Q c ) mo d e ■ 
While the difference between apparent and true V sin i val¬ 
ues seems to be significant in some cases, the increase of 
f l/H r is n ot so h igh, as compared to the ratio obtained by 
IChauville et al. ( 120011 ) with apparent V sin i parameters. 
This is simply due to the fact that we need to consider 
also the increase from apparent to pnre values of the stellar 
fundamental parameters. These last identify objects that 
are more massive and less evolved than expected from the 
apparent parameters. Such increase also implies a signif¬ 
icant augmentation of the equatorial critical velocity 14 
which is barely overweighted by the respective true V sin i 
value. 


Table 8. Iteration of S2/f2 c and (fl/fl c ) m ode 


6.2.1. Comments on the solutions obtained 

We remind that the values of parameters displayed in 
Table El are ’modes’ of the respective distributions of 9 3 
possible solutions. They can be considered as representing 
the most probable configuration of fundamental parame¬ 
ters for the rotating star, as they are perceived from the 
entry set of apparent parameters with their uncertainties, 
seen through the models calculated with fastrot used 
to interpret them. They look like a sort of independent 
’’determinations” which also seem sometimes to under¬ 
estimate the effects due to fast rotation. Nevertheless, if 
we used only the apparent entry parameters without their 
uncertainties, the unique translation into pnre quantities 
that comes out would lead in many cases to AT e g, A log g 
and Ay sin i differences between pnre and apparent values 
that are higher than those issued from the values in Table 
El We calculated the f2/fi c and (fi/f2 c )mode using these 
simplified estimates of pnre parameters and true ysini, 
as well as the same iteration procedure for f2/f2 c . We ob¬ 
taine d thus , = 0.83, near the value put forward in 

IZorec et al. (2003) and (12/S2 c )mode = 0.88. This small in¬ 
crease of the 12 /fh . as compared to the value obtained by 
IChauville et al. 1200ill , is entirely ascribed to the higher 
V c values implied by the apparent —> pnre transformation 
of fundamental parameters. 


6.2.2. pnre and actual fundamental parameters 


In this paper, we have studied the relation between 
the apparent parameters and those that should repre¬ 
sent the actual stellar mass and its evolutionary stage. 
The translation obtained into pnre parameters was per¬ 
formed with the help of evolu tionary tracks calculated for 
non rotating stars ( Schaller et all . 1992 1 However, evo- 
lutiona ry tracks of rotating o bjects are somewhat dif - 
ferent ( Heerer fc Langed 1200(1 iMevnet k Maedeil EoO(3l. 
Depending on the case, they may produce slightly differ¬ 
ent mass estimates. These estimates are those which we 
might finally consider as a better approach to the actual 
stellar mass. Nevertheless, the mentioned differences are 
altogether small and do not modify whatsoever the results 
obtained in the preceding sections. Their ins and outs are 
out of the scope of the present work and will be widely 
discussed in another contribution. 


modes 


fi/Qc 

V/Vc 

f2/f2c 

V/Vc 

n / n c 

M2 

M3 

M4 

0.800 

0.708 

0.842 

0.756 

0.879 

0.507 

0.369 

0.271 

0.842 

0.712 

0.846 

0.753 

0.878 

0.513 

0.375 

0.277 

0.846 

0.713 

0.847 

0.746 

0.872 

0.515 

0.376 

0.277 

0.847 

0.713 

0.847 

0.751 

0.876 

0.515 

0.376 

0.278 

1.000 

0.782 

0.898 

0.726 

0.854 

0.627 

0.513 

0.428 

0.898 

0.725 

0.856 

0.673 

0.814 

0.533 

0.397 

0.300 

0.856 

0.715 

0.848 

0.733 

0.862 

0.518 

0.379 

0.281 

0.848 

0.713 

0.846 

0.732 

0.861 

0.514 

0.376 

0.278 

0.846 

0.713 

0.847 

0.751 

0.876 

0.514 

0.375 

0.277 


7. Conclusions 

In this paper we present the calculation code fastrot 
to calculate the effects carried by fast rotation on the fun¬ 
damental parameters of early-type stars. Particular atten¬ 
tion is payed to represent hydrogen Balmer, Hei and Mgn 
lines in the blue spectral region, which are currently used 
to infer from spectroscopic data the fundamental param¬ 
eters of early-type stars. We also calculated a number of 
C, N and O lines in the visible spectral range, which enter 
the CNO abundance determination. 
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We have calculated synthetic spectra with classical 
plane-parallel non-LTE model atmospheres using appar¬ 
ent fundamental parameters derived from H 7 , He 14471 
and Mgil4481 lines and compared their predictions for 
other spectral regions of rotating stars. We noted that the 
above apparent fundamental reliably reproduce also the 
remaining Balrner and He 1 lines in the visual. However, 
discrepant fits are expected for CNO lines. The noted dis¬ 
crepancies may thus lead to systematic overestimates of 
the N abundances, while the C abundances can be under¬ 
estimated. Regarding the O abundance, it depends on the 
lines used. They are sensitive to local formation conditions 
and reflect preferences either on the polar or equatorial 
formation regions. O abundances can then be under- or 
overestimated, depending on the lines used, stellar appar¬ 
ent configuration, if rotational effects are neglected. 

The fastest rotation in the main sequence are proba¬ 
bly Be stars. The strong gravitational darkening can lead 
then to sensitive Vsini underestimations, which in turn 
may imply that the hitherto undercritical rotation of these 
stars can be argued. We payed much attention then on 
the determination of the fundamental parameters of fast 
rotators, in particular on Vsini. We thus obtained that 
although the classical Vsini parameter do underestimate 
the stellar equatorial velocity of fast rotators, the fact that 
the pnrc stellar mass is higher and the evolutionary stage 
lower than expected from apparent quantities, causes that 
in most cases the increase in the estimate of Vsini does 
not overweights considerably the new V c . The study of a 
sample of 130 Be stars leads thus to a most probable ratio 
V e /V c ~ 0.75, or 0/fl c ~ 0.88, which still imply average 
undercritical surface rotation. 
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Fig. 9. Comparison between spectra computed using pnrc (full line) and apparent parameters (crosses). Pnrc spectra 
were computed for T° s = 20000 K, log<? 0 = 4.0 and f2/f2 c = 0.99. Three different inclination angles were considered. 
The values of the apparent parameters for each case are given in Tabled 
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Fig. 10. Comparison between spectra computed using pure (full line) and apparent parameters (crosses). Pnrc spectra 
were computed for T° s = 20000 K, logg 0 = 4.0 and fl/Q c = 0.99. Three different inclination angles were considered. 
Apparent values for each case are given in Table 0 













Y. Fremat et al.: Effects of fast rotation 


19 


Table 9: Apparent and pnrc fundamental parameters of the studied 
stars. An electronic version of the table will be made available at 
the CDS. 


Apparent parameters Mode of the pnrc parameters 


HD 

T e s ± e 

log g ± e 

Hsin fie 

T e fi ± la 

log g ± la 

Fsiniiler 

Vo ±1(7 

i°±la 

144 

11956±324 

3.383i0.055 

125±07 

12206A414 

3.403i0.067 

132±08 

309±20 

33.9±1.2 

4180 

14438±272 

3.284±0.041 

195±10 

15584A426 

3.399i0.062 

208±13 

332±21 

55.7±3.2 

5394 

26431±618 

3.800i0.070 

432±28 

30001A624 

3.990i0.065 

441±27 

577±36 

76.4±2.7 

6811 

12478±444 

3.092A0.068 

85±05 

12578±481 

3.110A0.073 

93±06 

274±19 

26.5±1.0 

10516 

25556A659 

3.899i0.078 

440±30 

29031A751 

4.168A0.087 

462±33 

590±42 

72.5±3.9 

11606 

19484A297 

3.653±0.037 

280±14 

21422A439 

3.774±0.055 

286±16 

446±26 

58.2±3.5 

18552 

12715±284 

3.456i0.048 

286±16 

14646A244 

3.776i0.041 

292±14 

369±19 

75.2±5.0 

19243 

21243A388 

3.270i0.043 

160±08 

22752A409 

3.297±0.043 

165±09 

387±21 

36.2±1.3 

20336 

18684A517 

3.865i0.072 

328±21 

20843A558 

4.050±0.076 

341±23 

487±32 

66.7±4.0 

22192 

15767A509 

3.465±0.074 

275±19 

17690A299 

3.763±0.042 

295±15 

397±20 

75.0±3.1 

22780 

13299±284 

3.381i0.046 

285±15 

15686A443 

3.572i0.066 

297±18 

366±23 

90.Oil.0 

23016 

12230±241 

3.922i0.047 

235±12 

13420A249 

4.012A0.047 

242±12 

410±21 

48.7A2.6 

23302 

12754A504 

3.368±0.078 

170±12 

13484A293 

3.412i0.047 

181±10 

320±18 

46.8±1.6 

23480 

13691±481 

3.629i0.076 

240±16 

14828A402 

3.792i0.064 

254±15 

391±24 

57.3±3.3 

23552 

12711±416 

3.764i0.070 

220±14 

13708A239 

3.807i0.043 

230±12 

387±20 

52.0A2.2 

23630 

12258±505 

3.047A0.076 

140±10 

12885A292 

3.095A0.045 

149±08 

274±15 

44.6±1.6 

23862 

12106A272 

3.937i0.052 

286±16 

13436A257 

4.197i0.050 

290±15 

420±22 

65.3A3.9 

24534 

25191A513 

3.618i0.056 

293±17 

26831A699 

3.854i0.081 

294±20 

480±34 

54.7±2.6 

25940 

16158±582 

3.572±0.084 

197±14 

17593A338 

3.789i0.049 

220±13 

386±21 

38.3A4.5 

28497 

26724A427 

4.200±0.046 

335±17 

28101A738 

4.222i0.088 

342±24 

631±45 

45.4A2.0 

30076 

20488±330 

3.772i0.041 

213±11 

21311A397 

3.792i0.049 

225±12 

456±25 

38.4A1.7 

32343 

16131±514 

3.798±0.077 

95±06 

16121A430 

3.790±0.065 

103±06 

410±26 

19.7A0.8 

33328 

21137A599 

3.447i0.075 

318±22 

24340A644 

3.608i0.076 

333±23 

440±30 

73.0A2.4 

35411 

23652A510 

3.711i0.059 

170±10 

25138A441 

3.719i0.046 

174±09 

484±26 

29.Oil.3 

35439 

22134A665 

3.920i0.087 

263±19 

23547A364 

3.948±0.041 

266±13 

513±26 

42.2i2.3 

36576 

22618A508 

3.804±0.062 

265±16 

24502A399 

3.839±0.043 

266±13 

487±25 

46.3i2.2 

37202 

19310±550 

3.732i0.075 

310±7 

21460A420 

3.910i0.060 

326±7 

466±13 

66.0i4.1 

37657 

17951A542 

3.659±0.076 

198±13 

19275A301 

3.686±0.039 

209±10 

427±21 

41.9il.9 

37795 

12963A203 

3.517i0.036 

180±09 

13695A437 

3.559i0.069 

192±12 

355±23 

44.9i2.1 

37967 

16543±264 

3.850i0.041 

210±10 

17558A589 

3.873i0.085 

228±16 

455±33 

42.1i2.1 

38010 

25826±448 

4.089A0.050 

370±20 

28015A555 

4.194A0.061 

377±22 

608±36 

53.0i3.4 

40978 

18227A571 

3.576i0.078 

200±14 

19605A550 

3.608i0.073 

211±14 

401±27 

43.3il.9 

41335 

20902A610 

3.886±0.081 

358±25 

23431A639 

4.076i0.082 

376±26 

520±36 

68.5i3.8 

42545 

15764A334 

3.849i0.053 

285±16 

17522A309 

3.976i0.046 

303±16 

460±24 

60.1i4.0 

44458 

22914A532 

3.600i0.063 

242±15 

24745A667 

3.640±0.079 

254±18 

449±31 

48.0i2.1 

45314 

29160A467 

3.949i0.044 

285±14 

31092A557 

3.968i0.053 

291±16 

594±33 

40.6il.9 

45542 

13711±526 

3.627i0.082 

217±15 

14651A229 

3.682i0.038 

226±11 

369±18 

52.5i2.6 

45725 

17810±455 

3.895±0.066 

330±20 

19601A459 

4.168A0.065 

345±21 

486±30 

66.6i3.9 

45910 

20659±484 

2.739i0.052 

240±14 

22983A471 

2.992i0.048 

254±14 

342±20 

76.5i3.5 

45995 

21198±632 

3.769i0.083 

255±18 

22674A423 

3.806±0.050 

260±14 

470±26 

46.8i2.1 

47054 

12336±481 

3.538±0.078 

219±15 

13292A355 

3.607±0.058 

227±13 

351±21 

57.2i2.2 

48917 

20341±435 

3.395i0.053 

205±12 

21348A569 

3.438i0.070 

212±14 

389±26 

45.4i2.0 

50013 

24627A590 

4.017A0.072 

243±15 

25790A713 

4.026A0.087 

244±17 

535±39 

37.3il.9 

50083 

20908A604 

3.678±0.078 

170±12 

21216A625 

3.689±0.081 

176±12 

452±32 

32.5il.4 

53974 

27808±759 

3.650i0.085 

100±17 

27624A577 

3.649i0.060 

107±06 

490±30 

16.5i0.6 

54309 

20859A397 

3.587±0.048 

195±10 

22235A543 

3.610±0.066 

201±13 

425±27 

38.4il.5 

56014 

18544A614 

3.287i0.081 

280±20 

21061A453 

3.514i0.055 

290±17 

389±23 

72.6i4.2 

56139 

19537A331 

3.615i0.042 

85±04 

19465A300 

3.612i0.037 

89±04 

403±20 

16.6i0.7 

57219 

19865A575 

3.874i0.078 

80±05 

19759A486 

3.862i0.066 

84±05 

443±28 

14.3i0.7 

58050 

19961A465 

3.934±0.063 

130±08 

19998A590 

3.926i0.081 

136±09 

481±34 

21.9il.2 

58343 

16531±409 

3.619i0.059 

43±02 

16389A437 

3.611i0.063 

49±03 

377±24 

9.7i0.3 

58715 

11772±344 

3.811i0.062 

230±14 

12769A395 

3.858i0.068 

231±14 

380±24 

54.5i2.9 
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58978 

24445±476 

4.153±0.057 
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1 26622=4717 4.203±0.088 

375±26 

615=444 

55.2±3.0 

60606 

18030±332 

3.662±0.045 

273=414 

19600=4414 

3.783=40.055 

285±16 

441±25 

58.2±3.6 

60848 

27685±671 

4.051±0.077 

247±16 

28491=4630 

4.057=40.070 

256±16 

582=437 

35.4=41.8 

63462 

26558±635 

3.598±0.070 

440=429 

30685=4728 

3.835=40.077 

514=431 

545=437 

90.Oil.0 

65875 

20205±532 

3.845±0.071 

153=410 

20358=4590 

3.844=40.079 

163±11 

447=431 

27.8±1.2 

68980 

25126±642 

3.951±0.077 

145=409 

25172=4492 

3.945=40.056 

152±08 

534=431 

22.2±1.2 

75311 

16943±429 

3.846±0.064 

268=416 

18513=4527 

4.089=40.077 

283±19 

467=431 

54.2±4.0 

77320 

19555i488 

3.871±0.067 

345=422 

21539=4637 

4.015=40.085 

349±25 

501±36 

68.7±3.6 

83953 

15187±329 

3.661±0.051 

260=414 

16616=4273 

3.788=40.041 

276±14 

407=421 

61.0±4.0 

86612 

16374±577 

3.580±0.083 

185=413 

17392=4381 

3.790=40.055 

206±12 

381=422 

43.3±4.2 
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